Aging, phase ordering and conformal invariance 
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In a variety of systems which exhibit aging, the two-time response function scales as R(t, s) « 
s _1_a f(t/s). We argue that dynamical scaling can be extended towards conformal invariance, 
obtaining thus the explicit form of the scaling function /. This quantitative prediction is confirmed 
in several spin systems, both for T < T c (phase ordering) and T = T c (non-equilibrium critical 
dynamics). The 2D and 3D Ising models with Glauber dynamics are studied numerically, while 
exact results are available for the spherical model with a non-conserved order parameter, both for 
short-ranged and long-ranged interactions, as well as for the mean-field spherical spin glass. 
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Aging phenomena are observed in a broad variety of 
systems with slow relaxation dynamics |l|. Aging be- 
haviour is known to be more fully revealed by two-time 
quantities, rather than by one-time quantities (see [|2||3|] 
for reviews) . The most commonly studied two-time quan- 
tities are the correlation function C(t, s) = (<f)(t)(j)(s)) 
and the response function R(t, s) — 5(cj)(t))/Sh(s), where 
</> is the order-parameter field at some point, and h is the 
conjugate local magnetic field, while s is the waiting time 
and t the observation time. 

Consider for dcfinitcncss, instead of a genuine glassy 
system, the situation of a ferromagnetic model, evolving 
at fixed temperature T from a disordered initial state. In 
the high-temperature phase (T > T c ), the relaxation time 
is small, so that the system relaxes rapidly to equilibrium, 
where C(t) and R(t) are time-translation invariant (they 
depend only on the difference r = t — s) and obey the 
fluctuation-dissipation theorem: TR{t) — — dC(r)/dT. 
In the low-temperature phase (T < T c ), where coarsening 
takes place, both C(t, s) and R(t, s) depend non-trivially 
on the ratio x = t/s in the self-similar regime of phase 
ordering Q] . Similar features are also observed in the late 
stages of critical dynamics (T = T c ) |§-§|. The distance 
from equilibrium of an aging system is usually character- 
ized by the fluctuation-dissipation ratio X{t, s) |9J,[lO|,p|,p| , 
such that T R(t, s) = X(t, s) dC(t, s)/ds. 

In the scaling regime where s and r = t—s are simulta- 
neously much larger than the microscopic time scale (set 
to unity), the scaling laws 



C(t,s) 



b F(t/s), R(t,s) 



- a f(t/s) (1) 



are found to hold for a broad range of models |^-Q . More- 
over, for 1 > 1, i.e., 1 < s < (, both scaling functions 
usually fall off as 



F(x) ~ f(x) ~ x- 



(2) 



where z is the dynamic critical exponent and A is the au- 
tocorrelation exponent uM- For a ferromagnetic model, 



with scalar non-conserved order parameter, (|l|) holds at 
criticality (T = T c ) §,§], with a = b = 2(5/ {uz) = 
(d— 2 + if) 1 ' z, where /3, v, r\ are static exponents, while z 
is known from equilibrium critical dynamics. In the 
phase-ordering regime (T < T c ), (Q) only holds in the 
aging regime, where C(t, s) decays from its plateau value 
<Zea = Afg_ to zero (M eq is the spontaneous magnetiza- 
tion). One has 6 = and F(l) = qea- There seems to 
be no general result for the response exponent: a = 1/2 
for the Glauber-Ising model, both in one dimension |i]|7]] 
and in higher dimensions |12[, while a — d/2 — 1 for the 
spherical model in dimension d > 2 |l^ , |l4| ,|s| . We stress 
that the specific scaling forms (|l],||) only hold for a com- 
pletely disordered initial state. If there are correlations in 
the initial state, similar but different scaling forms hold 

At present, there is no general principle to predict 
the form of correlation and response functions in non- 
equilibrium systems. On the other hand, for an equilib- 
rium critical point (where formally z = 1), scale invari- 
ance can be extended to conformal invariance [p7|-fi9|. 
Scale invariance implies that correlators transform co- 
variantly under dilatations, i.e., scale transformations 
which are spatially uniform. Conformal transformations 
are local scale transformations with a position-dependent 
dilatation factor b = b(r), but such that angles are con- 
served. This is already enough to fix the form of an equi- 
librium correlator at criticality in any space dimension. 
Furthermore, 2D conformal invariance yields exact values 
for the entire set of critical exponents, the exact form of 
all n-point correlators, a classification of the universality 
classes, and much more |2(j] (see |18|,[l9[ for reviews). 

Is a similar extension of the scale invariance of (|l|) 
also available for non-equilibrium systems? Indeed, we 
shall present evidence that this might be so. First, we 
argue [pTJ that the two-time response function R(t, s) 
should transform covariantly under the action of confor- 
mal transformations in time. This assumption is then 



shown to imply the scaling form 

R(t,s)^r (t/ S ) 1+a -^ z (t- 
i.e., /(a;)=roa; 1+a - A / z (a;-l)- 



(3) 



The scaling function f(x) is thus entirely fixed, up to the 
normalization constant ro, by the two exponents a and 
X/z entering (Q) and (||). This explicit scaling form (||) 
of the response function is the main result of this Letter. 
It is expected to hold throughout the aging regime, for 
any value of the ratio x = t/s > 1. The power law (0) is 
recovered for well-separated times (a; ^> 1) p^ |. 

The prediction (|3]) will be corroborated by numerical 
simulations in the 2D and 3D Ising model with Glau- 
ber dynamics. In addition, exact results for the spherical 
model with a non-conserved order parameter (including 
spatially long-ranged interactions and/or quenched disor- 
der) also reproduce (|) |J,0||[l|J2|||]. This confir- 
mation provides, for the first time, evidence for conformal 
invariance in noncquilibrium and aging phenomena. 

We now sketch the line of reasoning leading to (0). 
The full calculation will be given elsewhere pl[ . To 
begin, we ask: what space-time symmetries are consis- 
tent with dynamical scale invariance t — > b z t, r — * br, 
where z is the dynamical exponent pqjjlfl ? A similar 
question has already been successfully raised for equi- 
librium systems with strongly anisotropic scaling p5| . p(| . 
While in equilibrium systems the correlation functions (of 
quasiprimary operators |20[ ) are expected to transform in 
a simple way, for non-equilibrium systems it is rather the 
response functions which will take this role, as argued 
long ago p7J . We expect the requested extension of dy- 
namical scaling to contain the Mobius transformations of 
time: t >-> t' = (at + f3) / (jt + S) , with aS — /3-y — 1, since 
these occur in the two known special cases, namely con- 
formal invariance for z = 1 and Schrodinger invariance 
for z — 2. It turns out that this condition is already suf- 
ficient to fix the form of the infinitesimal generators. In 
one space dimension, to which we restrict for notational 
simplicity, it can be shown that for non-equilibrium sys- 
tems one may write |21| 



X X 



-t 2 d t 



X = -td t 

- (2/z)trd r 



(l/z)rd r , 
pr 2 d 2 ~ z , 



(4) 



where j3 is a constant related to "mass" |25|J21[| . The gen- 
erators X n satisfy the commutation relations [X ni X m ] = 
(n — ni)X n+m (with n, m G { — 1, 0, 1}) of the Lie algebra 
of the conformal group. 

Equation ([|) forms the basis for the derivation of (||). 
Time translations are generated by X-\. In order to 
make the above construction applicable to non-equilib- 
rium situations, we must discard the latter, and only re- 
quire covariance under the subalgebra S generated by X 
and X\ (see [E3| ) . It is clear from the form of the gen- 
erators that the initial line t — is invariant under the 



action of S. Now, consider a general response function 
G = ((j>i{t\,ri)4>(f2,r2)), where the field <\>\ is charac- 
terized by its scaling dimension x\ and "mass" /3i, and 
the response field 4>2 (see e.g. jl7|]) has scaling dimension 
X2 and "mass" (32- Then the covariance of G under the 
local scale transformations in S is expressed by the con- 
ditions |l) X G = (Ci + &G, X X G = (2&ti + 2( 2 t2)G, 
where Q — Xi/z (i — 1,2). Moreover, we require spatial 
translation invariance, thus G = G(ti,t 2 ;r 1 — r 2 ). This 
is always satisfied if (i\ + (— l) 2 ~ z /32 = 0. We can now 
set r — 7'i — ?*2 = and obtain the response function 
G = i?(7j 1 ,rj 2 )- Then the generators (0) reduce to the 
standard conformal generators (see e.g. [|18| — [20f| ) , and R 
satisfies the differential equations 



(t 2 d t 



(td t + sd s + ( 1 +{ 2 )R{t,s)=Q, 
+ s 2 d s + 2(^ + 2(28) R(t,s) =0, 



(5) 



hence R{t,s) = r (/J/s) C2_Cl (i - s)~ Cl " C2 . An identifica- 
tion of exponents with (|l|) , (g) yields (g) . 

The prediction (||) will now be checked against results 
for various model systems. We begin with a novel numer- 
ical investigation of the Ising model, on square or cubic 
lattices, with periodic boundary conditions, and Glauber 
or heat-bath dynamics. Because the instantaneous re- 
sponse function R{t, s) is too noisy to be measured in a 
simulation, we consider instead the integrated response 
function @| 

p{t, s)=T f du R(t, u) « (T/h) M TRM {t, s), (6) 
Jo 

where Mtrm(£, s) is the thermoremanent magnetization, 
i.e., the magnetization of the system at observation time t 
obtained after applying locally a small magnetic field h 
between the initial time t = and the waiting time t = s. 
This quantity can be readily measured, either in TRM 
experiments, or in numerical simulations |2S|]. The data 
shown in Figures yj and |2j have been obtained in this way, 
and averaged over at least 1000 different realizations of 
systems with 300 x 300 spins in 2D and 50 x 50 x 50 spins 
in 3D. Larger systems were also simulated, in order to 
check for finite-size effects. Table | contains the numeri- 
cal values of exponents used in the subsequent analysis. 
Figure fl] displays our results for p(t, s) in the scaling 
regime at criticality (data corresponding to 1 ~ r <C s 
have been discarded). From (|l|) we expect a data collapse 
if s a p(t, s) is plotted against x — t/s, and this is indeed 
the case. Having thus confirmed the expected scaling, we 
can compare with the prediction (H). We find complete 
quantitative agreement, after adjusting only one parame- 
ter, the overall normalization constant ro. Figure @ shows 
our results in 2D and 3D, in the scaling regime and for 
a fixed temperature below T c . As a = 1/2 |0|, we ex- 
pect that s 1 / 2 p(t, s) only depends on t/s. This is indeed 
the case in 3D, but for 2D the situation is more compli- 
cated. Analytical calculations in the spirit of the OJK 



approximation reveal the presence of extra logarithms: 
p(t,s) ~ s _1 / 2 lns f{t/s) iQ. Logarithmic corrections 
to scaling are not so rare, even in the realm of 2D equi- 
librium conformal theories |2{]|Pj] . We therefore propose 
the heuristic ansatz 



p{t, s) w s- 1 / 2 (r + rilns) E(t/s), 



(7) 



where tq, ri are non-universal constants, and E(x) is a 
scaling function. It is apparent from Figure @(a) that we 
thus obtain a satisfactory scaling. We have checked that 
the same scaling form holds in the entire low-temperature 
phase, where r$, r% depend on T. We again find complete 
agreement between the form of the scaling function E(x) 
and the prediction (|3|). 

We now turn to confirmations of (|3j) by means of avail- 
able analytical results. For the ferromagnetic spherical 
model ||], which can alternatively be described in terms 
of a continuum field theory |3,[l5| , the scaling expression 
of the response function has been derived in any dimen- 
sion d > 2. It reads R(t,s) ~ (Airs)~ d / 2 f{x), where the 
scaling function is, in the ordered phase (T < T c ) uBM, 



fix) 



- ~d/4 



(x-l)- d '\ 



and at the critical point (T = T c ) |||| 

m - 



x i-*/*( x _ l)-rf/2 (2<d<4), 
(a; - l)- d/2 (4 < d). 



(8) 



(9) 



These expressions are in full agreement with (ra). The 
second expression of (pi) , corresponding to the mean- field 
situation, coincides with the result for a free (Gaussian) 
field 113], as could be expected. 

The spherical model has the peculiarity that the dy- 
namical exponent z — 2 throughout. In that case, the 
response functions are expected to transform covariantly 
under the Schrodinger group | J23|| . The full space-time 
dependent response reads 



4>(t, ri)</>(s, r 2 ) ) = R(t, s) cxp - 



M(n 



1-2 



(10) 



with R(t,s) given by (||), and where the "mass" A4 is 
a constant (in (Q), (3 = M./2 for z = 2). A comparison 
with the exact spherical model results, both at and below 
T c HH, also permits to confirm this fully HHJ. 

Recently, correlation and response functions have been 
calculated exactly for the spherical model with long- 
range interactions of the form J(r) ~ |r|~ d ~ CT p3. For 
d > 2 and a > 2, the spherical model with short-ranged 
interactions, discussed above, is recovered. On the other 
hand, for d > 2 and < er < 2, or d < 2 and < a < d, 
the dynamical exponent reads z = a below criticality, 
while the response function scales as p4]| 



R(t,s)^r (t/ S ) d /^(t~s)- d ^, 



(11) 



which again agrees with (Q). This example illustrates 
that spatially long-ranged interactions need not destroy 
conformal invariance in non-equilibrium situations, in 
contrast to the situation of conformal invariance at 
equilibrium. Moreover, for the mean-field spherical 
spin glass [^3|, the response function reads R(t,s) ~ 
(t/s) 3 / 4 (t — s)~ 3 / 2 in the low-temperature aging regime. 
This result also agrees with (||). It coincides with (||) 
for d — 3, as a consequence of the known similarity be- 
tween the 31? spherical fcrromagnet and mean-field spin 
glass [Ol. Finally, note that for the simple random walk 
R(t, s) = r = cste. for t > s, see pj, which is consistent 
with (|3|) with exponents a = — 1 and X/z = 0. 

In conclusion, the dynamical scale invariance realized 
in non-equilibrium aging phenomena apparently general- 
izes towards (a subgroup of) conformal invariance. As a 
first consequence, we obtained the explicit scaling expres- 
sion (||) for the two-time response function R(t, s), whose 
functional form only depends on the values of the expo- 
nents a and X/z. Our prediction (|3|) has been checked 
against analytical and numerical results for several spin 
systems with a non-conserved order parameter. The en- 
tire evidence available at present comes from classical 
systems with a disordered initial state. Different ini- 
tial conditions may lead to a modified scaling behaviour 
]l5| , [l6| and the applicability of conformal invariance to 
these remains to be studied. The problem of identifying 
the full set of physical conditions on the systems which 
obey (pi) remains open. However, validity of (0) might 
extend to a broader class of systems than studied here, 
possibly including some realistic glassy systems. Finally, 
it appears that the conditions for the applicability of con- 
formal invariance in non-equilibrium situations (such as 
the value of z, or the presence of long-ranged interac- 
tions) are less restrictive than for conformal invariance 
at equilibrium critical points. We hope that the ideas 
presented might shed some light on some of the standing 
questions of aging phenomena |50|] . 
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FIG. 1. Scaling of the integrated response function p for 
(a) the 2D and (b) the 3D Glauber-Ising model at criticality 
(T = T c ). The symbols correspond to different waiting times. 
The full curve is the conformal invariance prediction for p, 
obtained by integrating (J3J). 
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FIG. 2. Scaling of the integrated response function p for the 
low-temperature Glauber-Ising model (a) in 2D at T = 1.5 
and (b) in 3D at T = 3. The symbols correspond to different 
waiting times. The full curve is obtained by integrating (m). 



TABLE I. Critical temperature and exponents of the 2D 
and 3D Glauber-Ising model, both in the phase-ordering 
regime (T < T c ) and for critical dynamics (T = T c ). 





2D 


3D 


T c 


2.2692 


4.5115 


z 


T = 
T = T C 


2 
2.17 


2 
2.04 


A 


T = 
T = T C 


1.25 
1.59 


1.50 

2.78 


p/v 


T = T C 


1/8 


0.517 



